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The notion of a two-point susceptibility kernel used to describe linear electromagnetic 
responses of dispersive continuous media in non-relativistic phenomena is generalized 
to accommodate the constraints required of a causal formulation in spacetimes with 
background gravitational fields. In particular the concepts of spatial material in- 
homogeneity and temporal non-stationarity are formulated within a fully covariant 
spacetime framework. This framework is illustrated by re-casting the Maxwell- Vlasov 
equations for a coUisionless plasma in a form that exposes a 2-point electromagnetic 
susceptibility kernel in spacetime. This permits the establishment of a perturbative 
scheme for non-stationary inhomogeneous plasma configurations. Explicit formulae 
for the perturbed kernel are derived in both the presence and absence of gravita- 
tion using the general solution to the relativistic equations of motion of the plasma 
constituents. In the absence of gravitation this permits an analysis of coUisionless 
damping in terms of a system of integral equations that reduce to standard Landau 
damping of Langmuir modes when the perturbation refers to a homogeneous station- 
ary plasma configuration. It is concluded that constitutive modelling in terms of a 
2-point susceptibility kernel in a covariant spacetime framework offers a natural ex- 
tension of standard non-relativistic descriptions of simple media and that its use for 
describing linear responses of more general dispersive media has wide applicability in 
relativistic plasma modelling. 
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I. INTRODUCTION 



The behaviour of a material medium in response to electromagnetic and gravitational 
fields encompasses a vast range of classical and quantum physics. For media composed 
of a large collection of molecular or ionized structures recourse to a statistical description 
is required and this often leads to a coarser description in terms of a few thermodynamic 
variables and their correlations. Such a description relies on the efficacy of particular consti- 
tutive models or phenomenological constitutive data that serve to circumscribe its domain 
of applicability. 

For phenomena where the relative motions of the constituents approach the speed of 
light in vacuo or the material experiences bulk accelerations or gravitational interactions 
such constitutive descriptions must be formulated within a relativistic framework. However 
even within a spacetime covariant formulation there remains great freedom in how to ac- 
commodate electromagnetic responses that depend on material dispersion induced by spatial 
correlations or temporal delays of electromagnetic interaction^. The incorporation of such 
effects in a theoretical description often relies on a detailed structural model of the medium 
particularly if it is inhomogeneous or external gravitational gradients are relevant. Notwith- 
standing these complexities simple constitutive models have proved of considerable value 
for homogeneous polarizable media that exhibit temporal dispersion in a laboratory frame 
where gravity plays no essential role. Indeed the notion of permittivity and permeability 
tensors is often adequate to parametrize a large range of experimental linear responses of 
simple polarizable media to external static and dynamic electromagnetic fields. More gen- 
erally, for non-dispersive media these tensors can be subsumed into a susceptibility kernel 
that readily accommodates special relativistic effects on the bulk motion of media. 

In this article the degree to which the notion of a susceptibility kernel can be generalized 
to describe linear electromagnetic responses of dispersive continuous media is explored. In 
particular the effects of spatial material inhomogeneity and non-stationarity will be formu- 
lated within a fully covariant spacetime framework. In this manner the formulation can 
accommodate arbitrary gravitational and electromagnetic interactions. The framework will 
be illustrated by re-casting the Maxwell- Vlasov equations for a collisionless plasma in a form 
that exposes a 2-pointP' electromagnetic susceptibility kernel in an arbitrary external grav- 
itational field. This permits the establishment of a perturbative scheme for non-stationary 
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inhomogeneous plasma configurations in terms of sucli a kernel. Explicit formulae for the 
perturbed kernel are derived in both the presence and absence of gravitation in terms of 
the general solution to the equations of motion of the plasma constituents. In the absence 
of gravitation this permits an analysis of collisionless damping in terms of a system of in- 
tegral equations that reduce to standard Landau damping of Langmuir modes when the 
perturbation refers to a homogeneous stationary plasma configuration. 

It is concluded that constitutive modelling in terms of a 2-point susceptibility kernel 
in a covariant spacetime framework offers a natural extension of standard non-relativistic 
descriptions of simple media and that its use for describing linear responses of more general 
dispersive media has wide applicability in relativistic plasma modelling. 

II. CONSTITUTIVE RELATIONS 

In the following spacetime M is considered a globally hyperbolic, topologically trivial 
four dimensional manifold endowed with a metric tensor g with signature (—1, +1, +1, +1) 
describing gravitation. A closed 2-form F describes the electromagnetic field. The bundle 
of exterior p— forms over M is denoted hPM and its sections FA^M are p— forms on M . The 
bundle of all forms is AM = |Jp~gA^M. Associated with g is the Hodge map -k. Thus for 
a G FA^M its corresponding Hodge dual is denoted -ka G F^~^AM. The tangent bundle 
over M is denoted TM and its sections TTM are vector fields on M. We call the 1-form 
J = g{J, — ) G FA^M the metric dual of the vector field J G TTM. Maxwell's equations for 
the electromagnetic field F G FA^M in a polarizable medium containing an electric current 
J G TTM, satisfying the continuity (or current conservation) equation d'kj = 0, are written 

dF = and d*G = - (1) 

The excitation 2-form G G FA^M can always be expressed 

G = eoF + U (2) 

in terms of the permittivity eo of free space. The polarization^ 2-form H G FA^M results 
from all electromagnetic field sources not made explicit in J. 

In general H and J are non-linear functionals of F and other fields such as matter and ini- 
tial data on any initial spacelike hypersurface Sm C M. Such functionals are the constitutive 
relations describing G and J in terms of F and these other fields. 
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It is convenient to introduce integration on a fibred manifold M of dimension n + r with 
projection 7r_^ : M ^ N over a manifold of dimension n. Thus at each point a E N 
one has the fibre Mo- = T^J/^i^^} = {{(^' ■, ^) ^ | t^n{(^\ ■?) = ^} so dim(A/'o-) = r is the fibre 
dimension. For a G FA^^'^A/' we defind^ the form i a G FA'^'A^ by 



/3 A i a = I 71^(13) A a (3) 

N Jttj^ Jm 

for all /3 G FA"-PAr. 

In terms of local coordinates (a^, . . . , a") and (a^, . . . , a"", . . . ^'') for patches on A^ and 
M respectively, one may write the fibre integral 



(la] = ^ da^^ A ... A da^" I 



'''d/dah ■ ■ ■ '''d/dah a I {(7,?) (4) 

where = ttJ^^{{o'}) is the fibre over the point a E N and iQiQ^'k is the contraction on 
forms. Observe that if a does not contain the factor dq^ A ■ ■ ■ A dq"^ then i a = 0. The 
proof of this is given in appendix lemma |2} 

A key result of fibre integration, used to establish the current continuity equation, is that 
it commutes with the exterior derivative: 



(i 4 a 



da 



(5) 



for cr not on the boundary of A^ provided the support of a does not intersect the boundary 
of M . The proof is given in appendix lemma [3] 

In general models for 11 demand a knowledge of the dynamics of sources responsible for 
polarization as well as any permanent polarization that may exist in the medium. A full 
dynamical description depends on a specification of appropriate initial value data <^ on Ejvf. 
The exact structure of <^ depends on the sources of the polarization. For the plasma model 



described in section |III| the initial data corresponds to the velocity profile for each particle 
species at each point on Em in the plasma. 

In this article 11 is considered to be an ajfine functional of F of the form 

n[F,c]= / xApUF) + z[C] (6) 

''PX 

for some functional Z of The first term on the right is expressed in terms of the fibre 
integral of a two-point susceptibility kernel x £ FA^(Mx x My) expressible locally as 

X = IXabcdix, y)dx'' A dx^ A dy" A dy'^ (7) 
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Here Mx and My are two copies of M , locally coordinated by . . . , x^) and (?/°, . . . , y^) 
respectively, with projections px '■ Mx x My — )■ Mx, py : Mx x My — )■ My, px{x, y) = x, 
Py{x, y) = y and initial hypersurfaces T^Mx ^ and C My. Throughout, summation 
is over Roman indices a,b,c = 0, 1, 2, 3 and Greek indices /i, i/, a = 1, 2, 3. 

To consistently remove any reference to M (without a subscript) let F G FA^My, eo-F G 
TA^Mx, G G TA^Mx, J G TTMx and n[F,C] e TA^Mx- Thus eo can be regarded as a 
map eo : FA^My — t- FA^Mx which is the puUback of the natural isomorphism Mx — )■ My, 
together with a scaling to accommodate the choice of electromagnetic units. 

In terms of local coordinate bases on Mx and My the components of (|6| are 

U[F, CUix) = [ Ixabcdix, y) F,f{y) dy"''^ + Z[CU (8) 
in a multi-index notation with 

and 

(Note the reverse order for internal contraction.) Summations over multi-indices / C 
{1, . . . ,n} considered as an ordered p-list Ii < I2 <■■■< Ip of length |/| = p will also 
be employed. Thus 

dx^ = dx^"^'"^" = dx^^ A ■ ■ ■ A dx^^ 

and 

■ {x) _ .{x) _ ■ 

SO that, via summation, if a G FA^M then dx^ A i^f^a = a where |/| = p. 

In this notation the product manifold Mx x My inherits the following maps that will be 
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employed below: 

dx : TA^iMx x My) ^ rAP+^(Mx x My) , 

(i^(a) = ^^^dx" A dx^ A du-^ 
ax" 

dy : rAP(Mx X My) ^ rAP+^(Mx x My) , 
dyia) = ^S^dy" A rfx^ A dy^ 

-kx ■■ TA{Mx X My) ^ rA(Mx x My) , 
*x (a) = aij{-kdx^) A 

where a = ajj dx^ A dy^ 

Since F = (iA and for A with compact support away from any boundary of My it follows 
from ^ that 

n[F,c] = -/ (rfyx)Ap^y(A) + z[c] 

Hence n[F, remains invarianiP under the gauge transformation 

X X + dvC (9) 

for any ( = (abcdx""^ A dy'^ G rA^(Afx x My). Since the support of A can be made arbitrarily 
small dyX is uniquely specified by 11 [F, <^]. Furthermore 

rf*n[F,C] = - / (dx^xdrx) Ap*Y{A) + d^Z[C] 

•Jpx 

hence d-kI\.[F, ^] is invariant under the gauge transformation 

X X + dvC + ^xdxi (10) 

for any Q = (abcdx"'^ Ady'^ and ^ = ^abcdx"' A dy^^ . Similarly dx^^xdyX is uniquely determined 
hYd^Yl[FX]- 

In general, the permittivity functional 11 is a non-local functional in spacetime given by 
the integral ([s]). If x is smooth, and not identically zero, then 11 is always non-local. However 
for distributional susceptibility kernels it is possible for H to remain local. In this category 
one has the local, linear Minkowski constitutive relations 

H[F] = eo(e, - l)i,F Av + eoifi;^ - 1) ^ {{i, ^ F) A F) 
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where v 6 TTMy is a vector field representing the bulk 4-velocity of the medium and 
e^, fir £ TA^My are the relative permittivity and permeability scalars of the medium. These 
relations can be represented by a distributional susceptibility kernel with support on the 
diagonal set {{x,y) G Mx x My|x = y}. 

In general 11 is said to be causal on all of M if 11 1^ only depends of the values of F 
which lie on or within the past light-con^^ J~{x) C My of x. If 11 depends on <^ it may 
be causal on where M-J = T,Mx U li^s to the future of T^Mx}- The functional 11 is 
causal on if 11 [F, <^]|^ only depends on the values of F and <^ which lie on or within its 
past light-cone J~{x) fl M'^ of x and x G M'^. The data functional Z is casual on if 
^[C]U depends only on <^ G Sa/x ^ J {^) ^ ^ ^x- For 11 to be causal on M'^ it is 

necessary and sufficient (lemma |5] in the appendix) that the following be satisfied: 

• Z is causal on Mj^, 

• {dyx)\{x,y) = for all {x,y) G x My such that y ^ J~{x) and 

• 4My(^) I i^^y) = for all {x,y) G M+ x Sa/^ such that y ^ where is^,^ : 

X TjMy ^ Mx X My is the natural embedding. 

A. Spacetime homogeneous constitutive relations for media in Minkowski 
spacetime 

Minkowski spacetime has properties that underpin the notions of material spatial homo- 
geneity and stationary processes. Being isomorphic to a real 4-dimensional vector space it 
can be given an affine structure in addition to its light-cone structure. Physically this im- 
plies that no particular point in a spacetime without gravitation has a distinguished status 
and the concepts of material and field energy, momentum and angular momentum can be 
defined in terms of the Killing symmetries of the spacetime metric. Since all points of the 
spacetime are equivalent relative to this affine structure it is sufficient to denote Mx and 
My by M and, relative to any point chosen as origin, a point with coordinates x can be 
identified with a vector denoted by s G M^. It is then convenient to introduce the Minkowski 
translation map : M — )■ M, Az{x) = x + z that maps points a; to a; + 2; on M. 

If the electromagnetic properties of an unbounded medium are independent of location 
in spacetime they will be called spacetime homogeneous. Such electromagnetic constitutive 
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properties imply that variations in F at event y & M produce an induced variation in a func- 
tional nni-F] at event x G M, via a kernel Xabcd{x, y) that depends on the 4- vector x — y. If 
the constitutive relation is causal then there is no induced variation if x ^ J^{y). Further- 
more in a spacetime homogeneous medium Z[{^] = Zu where Zu G FA^M is independent of 

C- 

In terms of an electromagnetic constitutive functional IIh is given by 

Ilu[F]= I xAPrin + Zn (11) 
The functional IIh is said to be spacetime homogeneou^ if 

U^iA^F] = A^UulF] (12) 
This follows if the susceptibility kernel x satisfies 

X\{x+z,y+z) = X\{x,y) (13) 

and A*Zii = Zn. The contribution Z^ may model the presence of an externally prescribed 



stationary uniform permanent magnetic or electric polarization. Equation (13) implies the 
components of x in ([T]) can be written 

Xabcd{x, y) = Xabcd{x - y) (14) 

where 

Xabcd{x) = Xabcd{x, 0) (15) 

Thus, in a Minkowski spacetime for materials with electromagnetic spacetime homogeneous 
properties, ^ can be written in terms of a convolution integral: 

lln[FU{x) = \ [ X,Ux - y)F,f{y)dy''''f + {Z^),, 

JyeM (16) 

= le'^^'^iXabcd* F,f){x) + {Zn)ab 

where e'^'^^^ = ±1,0 denotes the Levi-Civita alternating symbol in coordinates in which the 
metric tensor takes the form g = rjabdx"' ® dx^ where rjab = diag(— 1, +1, +1, +1). In these 
coordinates the {Ziijab are all constants. 
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Let Fefik) and ^■H[F]ab{k) denote the Fourier transforms of Fef{x) and nH[-F]ab(a;) re- 
spectively, i.e. 



Fefik) = [ Fef{x)e"'-^dx''^' 



and 



where k = kadx"", k ■ x = kaX"'. Similarly let Xab^^ [k) be the Fourier transformation of 
U'^^fXabeaix), i.e. 



Xab'f{k) = [ XaU^)e'^-^dx'^^^ (17) 



If Zh = then it follows from (16) that: 

i^ii[FU{k) = \Xab'\k) Fe,{k) (18) 
Since Xabcd is a real function on M its Fourier transform satisfies 

Xab'^k)* = Xab'\-k) 

The 36 components of Xab'^'^{k) subject to this symmetry can be expressed in terms of 
permittivity, permeability and magneto-electric tensors relative to any observer frame. A 
specification of these components together with relations that determine the electric current 
J serve as an electromagnetic model for a spacetime homogeneous medium in Minkowski 
spacetime. If the medium lacks this electromagnetic homogeneity recourse to the Fourier 



transform (16) is not possible and the constitutive properties must be given in terms of a 



2-point kernel and ([8 



III. CONSTITUTIVE MODELS FOR A COLLISIONLESS IONIZED 
PLASMA 

As noted in the introduction the computation of the susceptibility for homogeneous sta- 
tionary dispersive media owes much to phenomenological models and input from experiment. 
For certain conductors, semi-conductors, insulators and low-dimensional structures much can 
also be learnt from the application of quantum theory. For inhomogeneous and anisotropic 
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media subject to non-stationary electromagnetic fields linear responses are often the subject 
of a perturbation approach. This is particularly so in the case of ionized gases. 

As an application of the above formalism the classical linear response of a fully ionized 
inhomogeneous non-stationary collisionless plasma to a perturbation is considered in the 
presence of an arbitrary background gravitational field. The perturbed constitutive tensor 
will be calculated in terms of solutions to the classical Maxwell- Vlasov equations for the 
system. This system is described in terms of the electromagnetic 2-form F G FA^M"*" 
over a gravitational spacetime M"*", lying in the future of an initial hypersurface S^^, and 
a collection of one-particle "distribution" forms (of degree 6), 9'^°^ G rA^£^+ (one for each 
charged species of particle [a] with mass m^"^ and charge q^'^^) on the upper unit hyperboloid 
bundle it : £^ — )■ over M+. The 7-dimensional manifold £^ is a sub-bundle of the 8- 
dimensional tangent bundle TM^ over M+ whose sections are all future pointing time-like 
unit vector fields on M"*". Thus generic elements of £^ can be written {z^w) with z G , 
7r(2;, w) = z and g{w, w) = —1. The initial values of the one-particle forms are given on the 
hypersurface S^- where S^- = n^^lTim} C 

The Maxwell- Vlasov system is usually written in terms of the Maxwell system in vacuo 
and all sources are contained in the total current J G TTM^. This in turn is given by the 
sum over each species current 

where J^"^ G FTM"^. Thus in terms of F and J the Maxwell subsystem is 

dF = and eodi<F = -i<J (20) 

The dynamic equations for each O^'^ can be written succinctly in terms of forms on £^ 
and a collection of Liouville vector fields W^'^ G TT£^ describing the flow of the charged 
particles associated with each species [a\. 

= H^,,^){z,w) + ^V(,,^)(2'J^) (21) 

in terms of certain horizontal and vertical lift^. With these vector fields the distribution 
forms 9^°^ are defined to satisfy the collisionless conditions: 

c/^L-i = (22) 
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and 



(23) 



To close this system one requires: 



^JN = J QV^^ (24) 

J n 

The closure of 6"°^ leads, from ([s]), to the continuity equation for each species current: 

di.J^"'^ =d( le^A= I de^"'^ = (25) 



so the total current 3-form -kj is closed away from the boundary Em- 

A local coordinate system {z^, . . . , z^) for a region containing z on M"*" induces a local 
coordinate system {z^, . . . ,z^,w^,w'^,w^) on S^. Since C TM'^ the tangent vector for a 
generic element {z, w) G may be written 

d 



{z, w) = w" 



dz'^ 



e ST c T,M^ 



where = n~^{{z}) is the 3-dimensional fibre of over z coordinated by {w^ , w"^ , w^) 
and w^{z,w) is the solution to gatw^w'' = —1 with > 0. All indices in the range 0, 1, 2, 3 
are raised and lowered using g""^ and gab so that Wq = w"'gao- Given a pair of vectors 
{z,w), {z,v) E £T C TzM^ the horizontal lift of the vector {z^v) to the point {z,w) G £^ 
will be denoted 'H(z,w){z,v) G T(^z,w)£^ and is given by 



{z,w) 



(26) 



where r°e/ are the Christoffel symbols determined by the metric components g"'^. Further- 
more if g{v, w) = then the vertical lift of the vector {z, v) to the point {z, w) G £^ is given 
by 

d 



{z,w) 



£ T{z,w)£ 



(27) 



Thus from (21), each Liouville vector field in these coordinates can be expressed as 

d 



\{z,w) 



d_ 

dz"- 'V - ejv~y- - ' ''J J dw 

Denote by f2 G TA^S^ the natural 7-form measure on given in these coordinates by 



+ ( - r\f{z)w'w^ + ^Fef{z)g'''w^) 
dz" V ^ mL"i J dw" 



(28) 



Wo 



(29) 
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In ref. [TTl eqn. (94) it is shown that for all species \a\ 

di^^id\fl = (30) 
The distribution function /^"^ G TA^S'^ relative to Q for the species \a\ is defined implicitly 



via 



From (30, 31) it follows that (23) is equivalent to 



(31) 



(32) 



and from (24) the components of the species current \a\ are given in terms of f^^^Zjw) by 

(33) 



= ^1(det^)(z)|V^ ^,,^^^^^^^,,3 



wo{z,w) 

A. Perturbation analysis 

Let e^f e TA^£+ and Fi e TA'^M+ be perturbations of ^^"^ and Fq, i.e. 

glc.^ ^ gi^] + + . . and F = Fq + Fi + . . . 



where 



0, 



and 



ciFo = 0, eoci*Fo = - / ^ w 



gN 



(34) 



(35) 



i.e. given by substituting F = Fq into (28). Substituting F into (21) yields 



+ . . . where VT/"^ = Wt\Fi) and the map Wi : rA2M+ ^ rr£+ is given by 



w^r(i^i)i(.,-) = ±^V(.,.)(^(.,^)Fi) 



m 



The first order linear system for the perturbation (6*1, Fi) is then 



de^f = , 



rfFi = 
eod ik Fi = — 



[q"" for 



(36) 



(37) 

(38) 
(39) 
(40) 
(41) 
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U! 



FIG. 1. A segment of the solution curve Ci^x,v) to the unperturbed Lorentz force equation (46) 
with final position x, final velocity {x,v), initial position y = C(2;^^)(r) and initial velocity {y,u) = 

Using ([5]) and (38) it follows that each species current in the sum on the right hand side 



of (41) is closed away from the initial hypersurface T,m- In terms of the excitation field 



Gi e TA'^M+ equation (41) will be written 



ci^Gi = 



(42) 



where 



Gi = eoFi + ni[Fi,Ci] 



for some linear functional Hi of Fi and <^ such that 



(43) 



(44) 



and Ci = {C}"'^ , Ci"^ , • • • } where C}''^ = i'f \^^^ for some i'f G TA^£+ which solves 9^^ = di^f . 
Thus Ci"^ is related to the initial velocity profile of the species \a\ . 



In the next section IIIB the general susceptibility kernel x ^ ThP{Ml^ x My) and linear 
functional Zi, determined by 6^^ and Fq, are found such that 



ni[Fi,Ci]U.= / xAp^(Fi) + Zi[Ci] 



(45) 



satisfies (44). 



B. A general formula for the functional Hi in an unbounded plasma 



In this section a general expression for a susceptibility kernel will be constructed in terms 
of the integral curves of the vector field G TTS^ . Such curves describe segments of 
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particle world lines under the influence of the Lorentz force due to the external electromag- 
netic field Fq. Although, for a general Fq, it is not possible to derive an analytic form for 
such integral curves, special cases are amenable to an analytic analysis. 

It proves convenient to let the final and initial states of each species of particle reside in 
fibres over and My respectively, bounded by the equivalent hypersurfaces T,Mx C 
and T^My C My. Thus the corresponding upper unit hyperboloid bundles nx '■ £x ^ 
and TTy : £^y — 7- My with boundary hypersurfaces S^-^ C Sj^ and S^-^ C Sy are used to 
accommodate the final and initial 4-velocities of the particles. The generic elements of these 
bundles are written G £x (z/; ^) ^ where x G Mj^, y G My and g{v,v) = 

g{u, u) = —1. The induced coordinate systems for Sx and Sy are . . . ,x^, v^, v"^, f^) and 
(?/°, . . . ,y^, u^,u'^, u^). Let f v), vo{x, v), u^{y, u) and uoiy^ u) be defined in the same way 
as w^{z, w) and wq{z, w). 

The contribution to the tensor ni[Fi,<^]^] due to all dynamic sources, arises from all 
particle histories in the past light cone of a; G M^- The history of the species particle [a] 
which passes through event x with 4- velocity (x, v) G £x will therefore be parametrized by 
negative proper time r: C|°^^) : [tq"^ (x, f ), 0] — )■ M+, r i— )• Such a history is the 

unique solution to the Lorentz force equation 




[a] 




) 



(46) 



with 



gicf = -1 



(47) 



and final condition 



{x,v) 



(0)=^, C^l^{Q) = {x,v) 



(48) 



where Cg,)(r) = C£.).(9.|.) = Cg^^) 



ax' 



and the value Tq (x, f ) < solves 



(49) 




The family of a// such histories is described in terms of the maps 



0W : ATjfl ^ £^ 



Y ' 



(50) 
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where 



= { (r, X, v) eR- xSj^l T^"'^ {x,v) <T<0} 

The manifold Af^ with boundary is naturally a fibre bundle over with projection zu^^ : 
Afx — ^ ^x, (t, X, v) (-7- zj^ (r, X, v) = (x, f ) and for any form a G TA^Afx it follows from Q 
that 



a*-"'"''(r, X, t>)(ir 



where a = a*^^^ (r, x, f )(ix^ A c??/'^ A dr + a'^^^ (r, x, f )dx^ A dy^ . 

Let rA|g be the set of sections over S^:^ with values in A^^y , i.e. if a G rA|^ 
then for each {y,u) G S^^, G A^^^^£^y. Let the map ^9^"^ : rA|^ Sy — )■ FA^f^^f be 

given by 



(51) 



where (j)]^'^ : Sx , (p;^'(x,v 



j?^ (x, = (/)(r, X, f ). For each species [a] let the initial data be given 
by e FA|^^^+ with z^^HCr^ = 0. 

In terms of these maps, it will now be shown that the general polarization functional Hi 
on Mx is given by 

1_T JlTX 



(52) 



where Hi and Zi are arbitrary linear functionals of Fi and C,i respectively. The excitation 
ni[Fi,<^^], in (52), is the general solution to (44) where the source di satisfies (|38||39|). The 



first two terms on the right hand side of (52) are linear functionals of Fi whereas the last 



term is a linear functional of the initial data C,i- Clearly 7kr(i(Si[Fi]) and T^r(i(Zi[<^]^]) are in 



the kernel of d-k, the homogeneous differential operator associated with (44). 



The proof that (52) solves (44) requires the following lemma which is proved in the 
appendix. 

Lemma 1. Let N he a manifold with a boundary G N and let V G FTA^ be a non- 
vanishing vector field on N such that every integral curve ofV intersects S^v precisely once. 
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For each a & N let the integral curve ofV terminating at a be given by 7o- : [ro(o"),0] -> 
where 7,^(0) = a and 7^(ro(cr)) G Sjy. The set M = { (cr, r) C x iV | rmin(o-) < t" < O} 
is a fibred manifold over N with projection wj^j : ^ N, (r, cr) i— )■ Q7Ar(r, a) = cr. The 
family of integral curves of V can be described by the map (pN '■ Af ^ N, (f)]\f{T,cr) = 7o-(r). 
Let ( G TA^^N such that iyC = 0, i.e. ( is a p-form on with values in hPN . Let 
ipN : Tkl^N ^ Vk^N be given by ^.^(OU = 'pNM^CUa)) G A^^iV. 

If P E rA^N is a p-form on N with compact support such that iyP = and G FA^A^ 
has the form 



then 



(53) 



iyd^ = (3 



(54) 



and 



c. 



This lemma is applied with = Sj^, wn 
= C = ([--^ and 



zu 



X ' 



V = W^''\ To = r, 



5 'PA' 



/5 - ^Tyfl(Fi)^0^ 



(55) 



Thus ^ in (53) becomes the 5-form ^|°^ G FA^^^f; 



(56) 



smce deg (0^"^*(^vkW(Fi)^o'^)) = 5- o^^^^ *o satisfy (|38| let 



(57) 



Furthermore from (54) and (55) 



SO (39) is satisfied. In terms of .^j"^ (52) can be written 

Ili[F,Xi]\. = J2l"'^ / +*rf(Si[Fi])U + ^d(Zi[Ci]) 
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Then from ([s]) 



Thus the Maxwell equation (44) is also satisfied. That (52) is the general solution to (44) 



follows from the fact that the difference between any two solutions of (44) satisfies the 



homogeneous differential equation associated with (44). 

Thus we have succeeded in eliminating 6'}"^ from the perturbation system (38 41 ), thereby 
reducing the system to dFi = and 







(58) 



in terms of (6*0, Fq), for the perturbation Fi. The perturbation 6i is then given by ([57p6|. 



C. The susceptibility kernel for an unbounded collisionless plasma 



Equating (52) and (45) with the initial data 



(59) 



yields 



/ X^Py{Fi) = J2^"'^ f /,-rfrA0L"l^(z^L.l(^^)O+^4Si[^i]) (60) 



X M+ U M+ X S 



Away from the initial hypersurface boundary d{M^ x My) = T^Mx 
using ^ and (A2) one has 

/ ^xdxiAp*yiFi) = I *xdiApUFi) = ^d I ^ A p^(Fi) = 

Jvx Jvy J Vx 



where is a linear functional of Fi. The gauge freedom x ~^ ^x^x^ given in (10) is 



equivalent to the addition of the term •kd{^i\Fi\) in (52). 



If Fi is restricted to have support in a certain domain one may find x such that 



ale 



(61) 
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To find such a susceptibility kernel requires the following maps. 

For {y,u) G let C^"^^^^ : M+ ^ M+ and : [0, rf^ (y, m)) ^ S+ be the unique 

solutions to the unperturbed Lorentz force equation ( 46p7 ) with initial conditions 



Citu)iO)=y and Cg,.)(0) = (2/,n) 



(62) 



where T{''{y,u) G M"*" U {c>o} is the supremum of the values of r such that C^^^{t) G M. 



Let <l>L"l : AT,^"^ ^ M+ x M+, 



(63) 



where 

Af^' = {{r,y,u) G M+ X ^+ I < r < rfi(y,M)} 

This map gives the final and initial positions of a solution to the unperturbed Lorentz 
force equation in terms of the initial position, velocity and proper time parameter r G 

Observe that is never surjective, since if ^^"\T,y,u) = {x,y) then x G J^{y). Also 
is never injective since (0, y, u) = {y, y) for all (y, u) E Sy- Thus does not possess 
an inverse and one must work locally on x My in order to establish the diffeomorphism 



(64) 



I.e. 



M/L'^i(C(,,„)(r),i/) = (r,2/,w) 



with P C M+ X M+ and C N^^ given by 



and 



There exists a unique u G and r G M+ such that N(r) = X for all [a] > 



(65) 



= I (r, m) $^"1 (r, y,u) eV for all [a] | 
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This map encodes the solution to the two-point problem, namely given an initial event 
y G My and final event x G Mx find the unique worldline to the unperturbed Lorentz force 
equation which passes though these two points. This worldline is specified by its initial 
velocity {y^u) G £x ^"^^ proper time r. The statement that does not have an inverse 
is equivalent to the statement that in general there may not be a unique solution to the two 
point problem on an arbitrary domain. The domain V is the set of all pairs (x, y) such that 
there is a unique worldline. 
Set 

where 

x'-^^y) = I'^^.^xdy'^' ^^^:l,^^-^\dT (67) 

for points {x,y) G T). In the appendix (lemma |6| it is shown that given x G (61) and 
Fi with support in 

V^ = Vr}px^{x} = {yeMY\{x,y)eV} (68) 



then (61) holds at x. Furthermore although {dYx)\(x,y) is unique, x has the gauge freedom 



given by ([9]). 

One may write (67) implicitly as 

Ap^y7 = -q"-' S^'-'*{dT A ^V"^^(V^i(,)0) (69) 
for all 7 G TK^M+ where S : ^•%^y){M+ x M+) -> A|',,j,)(M+ x M+), 

S{a) = 4^23" A dy^'^' (70) 

The tensor projector S has the simplest representation in the coordinate basis employed 
here since i'^a^'dy^ = 6'^. 



From (64) for a chosen species [a] one must consider r and u to be functions of {x,y) 
as well as the species label [a]. Thus let \E'^"^ be given by the functions r = T{x,y) and 

= u^{x,y), where we have dropped the species label, i.e. T{x,y) and u^{x,y) solve the 
implicit equation 

^SU..))(^(^'^)) = ^ (71) 
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where vP{x,y) is the solution to u"'{x,y)u^{x,y)gab{y) = —1 and uo{x,y) = gao{y)u"'{x,y). 
Let fif = flf\y, u) represent the unperturbed probabihty function on Ey. The contribution 
to the susceptibihty kernel from species \a\ is given in local coordinates by (lemma [T] in 
appendix.) 

/rLal2 I Jp-j- ^|3/2 
X \{x,y) - Jo g U € ecbihe^^uaX 

\ f72) 

2 (9?/" dx"^ dx'' 2 dx'^ dy" dx'' 2 dx'^ dx"" dy" ^ ' 

+ ( - r%,u^uf + ^Fo,,9'^'^f)——j dx,, A df'^ 

where g, Fq and F^e/ are all evaluated at y E My and each r and u belongs to the species 
\a\ . This is a key result of our article. 



D. A spacetime inhomogeneous microscopically neutral plasma. 



In a Vlasov model, a plasma or gas is deemed microscopically neutral if in its unperturbed 
state Fq = 0. Let M be Minkowski spacetime with global Lorentzian coordinates so that 



V^f^ = 0. Assume that /q solves the zeroth order Maxwell- Vlasov system (35) with 9q = 
i w if^^'^Q) and Fq = 0. In this scenario one can calculate x explicitly. 

Since Minkowski spacetime is flat and Fq = the integral curves C(^x,v) in global Lorentzian 
coordinates are the straight lines: 



r = \/—g{x — y,x ~ y) and u 



{x-y) 



T 



(73) 



Differentiating with respect to x"" and gives. 



dx° 



dr 



and 



Qya ' Qj.b 

dy^ T 



(74) 



If follows from ( 72 ) that 



^ \ (x,y) 



Lqt^ 9'"'''^cMh{2dxQ^ + e''''^%,^u''uadx,k) A df'^ 



(75) 



where T{x,y) and u{x,y) are given by (73). 
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It is often useful to explore the response of an inhomogeneous plasma due to a monochro- 
matic electromagnetic plane wave with constant amplitude E: 



(76) 



Setting the initial hypersurface as S^:^ = {?/° = i/q}, the general initial 5-form G rA|.^ S- 
satisfying iwoCi^ = is given in terms of its components by 



Y 



C["^\{o,y^,u-^) = {u^dy^ - u^dy^) A {C^^dy^ A du^^^ + C^^dy^ A rfw^^s^ ^ C^^dy^^ A du 



(77) 



where = Ci\{y^,u'^) for A = 1, . . . 6. For the integral curves (73) and the initial 
hypersurface S^-^ = {y^ = ^/q} one has tq{x,v) = (?/q — x^)/v^ and the map ip is given by 



(51) with 0*(?/") = + ry" and 0*(m") = v"- From (45) with x given by (75) and 2'i[<^;^] 



given by (59) one has: 

ni[Fi,Ci] = 



2 

- — Ee"''^^ 



+ dx^^ I dv'''''T^"^v 



Km 77n Ci> j+rfa;^^ / dv ^(Qj Ci> 



- CI' 5 



- 1 



+ ^rf(Hi[Fi])+^(i(Zi[Ci]) 
where / dv^'^^ denotes the triple integral operator JJJ^^ dv^"^^, t>*^ = -^/l + v^vt^, 

'(yO-xO)/^o 



(78) 



(79) 



and C = C(x^^;'^) = C}:!4C 
plane fronted. 

For the particular case of a plane fronted plasma distribution: 



/f°,t;^) in (78). This response is not in general 



flf\x,v) = h^\x',x\v')5{v')5{v') 



(80) 



21 



with initial data: 



(78) becomes the plane fronted 2-form 



ni[Fi,c]|. 

. „Lal 2 poo pO 

(81) 



describing the response of a spacetime inhomogeneous unbounded plasma to (76). 



E. Spacetime homogeneous unbounded plasmas 

The previous discussion simplifies considerably if the unperturbed plasmas is homoge- 
neous in space and time. In Minkowski spacetime M, an unbounded unperturbed plasma 
is deemed spacetime homogeneous if A*Fq = Fq and A*9jf = 9jf for all z E M where the 
translation map : M ^ M, Az{x) = x + z induces the map A 2^ S — y A^ — A^^. 
Such spacetime homogeneity implies that in all inertial frames the medium is stationary and 
spatially homogeneous in all directions. Such a spacetime homogeneous plasma will give rise 
to a spacetime homogeneous electromagnetic constitutive relation. In addition to the com- 
ponents (-Fo)afe with respect to an inertial frame being constant, the functions f^^^^x^v) are 
independent of event position x and can therefore be written /'^^^(f). 



In this scenario the Fourier transform (17) of the susceptibility kernel (18) for each species, 
is then given by 

J-oo J "^0 ^ / 

(82) 

where Fq is the 4x4 real matrix with components (-Fo)"?) = v"''^{Fo)cb generating the matrices 

^ ^ ' (83) 

XL«l^(r) = / D^-'Mt') dr' , L^-\%t) = g,,L^-''',{T)g'^ , 
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k-L^-'^v = kaL^"^''b{T)v'' and 



(84) 



The susceptibility kernel (82) can be shown to agree with the results of O'SuUivan and 
DerfleM 

Furthermore for a microscopically neutral spacetime homogeneous plasma with Fq = 0, 



Gi=0 and /ft"^(f) = h^;\v^)6{v'^)6{v^) it follows from (g and Q that for Im(w) > 



(85) 



The relativistic Landau damped dispersion relation for plane fronted Langmuir modes in an 
unperturbed spacetime homogeneous plasma arises by analytic continuation of the integral 



(85) to the lower-half complex u plane. 



F. Langmuir modes for an inhomogeneous unbounded plasma in Minkowski 
spacetime 

If the plasma is microscopically neutral but spacetime inhomogeneous in its unperturbed 



state the Landau dispersion relation corresponding to (85) becomes more involved. We 



define the generalized Langmuir sector to contain perturbations described by (81) but with 
the external polarization specified by set to zero. Since C}"^ = 0, ni[Fi,0] will be 



denoted ni[Fi]. Thus (43) with Gi = becomes 

eoFi = -ni[Fi] 



(86) 



Consider the case where planar inhomogeneities in a plasma composed of electrons and 
ions arise from the unperturbed spacetime inhomogeneous solution to the Maxwell- Vlasov 



system: (35 36) with Fq = and 



/^^ {x', x\ x\ x\ v\ v\ v') = fr x\x\ x\ v\ v\ v') 

= h(x'-'^,v'y{v')6{v^) 

where q^^ = —q^°^. 

For example one might consider 



-Lioril/^O ^1 2 ^3 ^,1 ^,2 ^,3\ 



(87) 
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where 74l-*°°l(x^) normalizes (87). Then f^°^ initially at x'' = represents a distribution of 



ions where, at each spatial point x^, the velocities belong to the 1-dimensional Maxwell- 



Jiittner distribution. In such a distribution the temperature 



ioril f^l 



and the number 



density of ions n^°^{x^) depend on position. It follows from (87) that f^^ also initially rep- 
resents a position dependent Maxwell- Jiittner distribution where and 



rp\ion\ (^jr^^ru^b^ / m^"^ . After the initial moment, the ions and electrons drift accord- 



ing to (87) and velocities do not remain in the Maxwell- Jiittner distributions. Alternatively 



(87) might describe a plasma composed of particles and anti-particles. 



In the theory of a spacetime homogeneous plasma u and k satisfy the transcendental 



dispersion relation (85). This relation contains an integral that is potentially singular. The 



Landau prescription circumvents this singularity by complexifying u and defining an analytic 
continuation for the integral in the complex u plane. 



Setting Hq^x'^ , , v^) = h[x^ — v^x^ jv^ ,v^^ in (80) yields (87) and (81) becomes 



ni[Fi]U 



dv^hix^ 



v^x^ 



(88) 



To compare with the results (85) given for the homogeneous case, consider the limit Uq — )■ 
— oo with Im(ci;) > 0. Furthermore for the non-evanescent modes considered here Im(/c) = 0. 



Thus (88) becomes 



Iii[F^]\x = -dx'''e^QiEe 



dv 



^h{x^ — v^x^/v^,v^) 



(89) 



where 



Je^2 



Jell 2 



+ 



In a spacetime inhomogeneous plasma there is no time-harmonic solution or associated 



transcendental dispersion relation between u and k. We therefore propose solving (86) with 
a longitudinal field Fi represented as the packet 



dx 



01 



dCj / dk E{u, k)e 



-iuix^+ikx^ 



(90) 
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Substituting (Isgl) and ([90| into ([86| yields 



du / dk E{oj, k)e 



-iCox +ikx^ 



— oo ^ — oo 



Ql dCj dk E{(2j, k)e 



dv 



dv 



Performing the inverse Fourier transform gives 
ATT^E{u,k) 

/OO POO POO POO 

/ dx' du dk E(c:;,A;)e'(-(^-^)^°+(^-'=)^') 
oo J —oo J —oo J —oo 

Since 

oo POO 

dx^ / rfxie^(-('^-'^)^°+(^-^)^')/i(x^ - v^xyv^,v^) = 2nh{k - k,v^)S{Ld - cu + v\k - k)/v^) 

-oo J — oo 

where 

/oo 
e-''''his,v^)ds 
-oo 

one has 



E{uj,k) 



27T 



dv' 



— hik - k, v^)6(uj -u + v\k- k)/v°) 
v^iuv^ + kv^y ^ ' 

(91) 



Since we restrict to non-evanescent modes k and k are real. For Eiuj^ k) to be non-zero 
one requires the argument of the 5-function to be zero. Since is real and therefore 
v^{k — k)/v^ is real it follows that although \m.{uj) > and Im(c2;) > the difference u — u 
is real. Furthermore from ui — u + v^{k — k)/v^ = it follows that |a; — < \k — k\. Thus 



(91 ) becomes 



where 



I{u, k, k) 



E{uj,k) = ^ / dk I{LO,k,k) 



doj E{Cj, k) — — t^h { k~k. 



(92) 



S{iJ,k,k) 



[Cok — kuy 



k-k 



{k-kY-{Cj-uY 



(93) 



and the contour of integration for uj in (93) is the straight line S{u,k, k) where lm{u 



Im(u;) > and —\k — k\ < Re(c2; — u) < \k — k\. Since {u — uoY < {k — k)"^ the arguments of 



h in (93) are always real and non-singular on S{u, k, k). 
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UJ - 



C-plane 



Im(a;) A 



uj—\k — k\* 



Im(a;) > 



UJ 



• 



uj + \k-k\ 



Re 



Wo 



uj — \k — k\t 



tUJ + \k — k\ 



FIG. 2. The upper contour denotes S{uj,k,k) when Im(a;) > for real k,k. The lower contour of 
integration is used when lm{uj) < for real k, k. 



L 


Im(a;) 

k\ 


UJ - C-plane 


branch cut 




j2 > 


Re(a;) 

o 



FIG. 3. Branch cuts in ui for I{uj, k, k). 



To accommodate the situation when E{uj, k) describes damped electromagnetic waves 



one must continue (93) to Im(a;) < for real k. However there is a double pole in the 



complex u plane at uj = = ku/k that coincides with S{uj,k, k) when Im(cj) = and 



\uj\ < \k\. To define an analytic continuation of (93) to Im(u;) < when |Re(ci;)| < \k\, we 
indent S{u, k, k) to encircle the pole in the standard manner and write the contour integral 
in terms of a principle part and associated residue, see figure |2} Such a continuation scheme 
gives rise to branches in the oj plane for J(a;, /c, k) as shown in figure [ij 
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This analytic continuation of (93) to Im(a;) < aquires the residue 



R{u, k, k) 



k\k\ duj\k 
k 



E 



(fc2 _ ^2)3/2 V k 



uk 



k ) h„n k — k 



Vk^ 



where hyi{k,v^) = ^{k,v^), Sk = k/\k\ and Sf^_j;. = {k — k)/\k — k\. In the case when 
\m.{u) = 0, the principle value of (93) is taken together with residue ^R{uj, k, k). Equation 



(92) then gives 
k 



E{uj, k) 
E{oj, k) 

E{u}, k) 



27r 
k 

2^ 



27r 



— oo 
oo 



I{u,k,k)dk if Im(ci;) > or |Re(a;)| > 

/oo 
R{u, k, k) dk 
-oo 

if Im(w) < and |Re(a;)| < \k\ 



VI{u,k,k) dk- 



ik 



R{u, k, k) dk 



if lm(u;) = and |Re(w)| < 

(94) 



where VI {u, k, k) in (94) refers to the principle part of (93) when Im(c<;) = and |Re(a;)| < 
\k\ and hence the pole at Cjq lies on the contour ^(ci;, k, k). Thus in each domain above, the 
perturbation E{u,k) must be determined by solving a non-standard integral equation. 



IV. CONCLUSIONS 



In this article a classical covariant description of electromagnetic interactions in contin- 
uous matter in an arbitrary background gravitational field has been formulated in terms 
of a polarization 2-form that enters into the macroscopic Maxwell equations. Linear dis- 
persive constitutive relations arise when this 2-form is expressed as an affine functional of 
the Maxwell 2-form with the aid of a 2-point susceptibility kernel. We have explored the 
constraints on this kernel imposed by causality requirements, spacetime Killing symmetries 
and local gauge freedoms. The formalism has been applied to an analysis of constitutive 
models for waves in coUisionless plasmas. In particular a formula for the linear susceptibility 
of a fully ionized inhomogeneous unbounded non-stationary coUisionless plasma to a pertur- 
bation in the presence of gravity has been given in terms of maps describing the dynamics 
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of the plasma. This formula has been elucidated by reference to both homogeneous and 
inhomogeneous perturbations in Minkowski spacetime. In the former case one recovers the 
standard Landau dispersion relation when perturbing Langmuir modes. In the latter case 
we have described a generalized damping mechanism for such modes that may arise when 
the unperturbed state is both inhomogeneous and non- stationary. Such a mechanism arises 
from the analytic continuation of an integral equation that replaces the Landau dispersion 
relation. 

It is concluded that the use of a covariant 2-point affine susceptibility kernel in describ- 
ing the electromagnetic response of dispersive media offers a modelling tool that naturally 
generalizes the use of permittivity and permeability tensors used to model electromagnetic 
interactions in non-rclativistic media. The formulation in terms of an arbitrary background 
spacetime metric offers potential applications in a number of astrophysical contexts involving 
electromagnetic fields in inhomogeneous or non-stationary plasmas 
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Appendix A: Proofs of results used used in the text. 



Lemma 2. Local representation of a in (4) from the implicit definition in equation (3) 



Proof. On a fibred manifold Af of dimension n + r with projection ttj^ : ^ N over a 
manifold of dimension n. Thus at each point a E N one has the fibre Af^- = '^J/{o'} = 
{(cr', <;) G A/" 1 7Tj\f{cr', <?) = cr] so dim(A/'o-) = r is the fibre dimension. Let (cr\ . . . , cr") and 
(cr^, . . . , cr", . . . ^'') be local coordinates for patches on A^ and Af respectively. 

Consider first the case when a G rAP"'"^A/' consists of a single component ai{a,<;)da^ A 
(i<j^ "^ with no sum on I. Hence explicit summation will be used in this particular proof. Set 
/ = {1, . . . , n}\I so that da^ A da' = ida^-" and let /3 = Ej then j3 A da' = 
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±(3jaida so that: 



T J(cr.c:)eAf 

y2 / Pj{a)da-^ A a/((T, <;)da^ A d<;^ 
/ f3j{(T)da-^ Ada^ Aai{a,q)dq^ 



A...r 

3j{a)da-^ A ai{a, q)da^ A rfc^-'" 
3j{a)do-^ A (icr^ A a;/((T, ^)cic^-'' 
l3j{a)da^ A da^ A a;/(o-, q)dq^-'' 
(3j{a)da^ Ada^ [ ai{a,<,)dq^-' 
V /" f3j{a)da^ Ada^ [ ai{a,<;)d<;^-' 

/ (3 A da^ I ai{a,<;)d(;^ 



A...r 



Thus by hnearity 

J JctGN JMa 



M 



where a = J2i '^/(c"; '^)da^ A dq^---'"' . If Q holds then for a = J2i '^/(c"; '^)da^ A dq^' 
Jn Jnj^ J Jn J<;eMa ' J Jn J^&Ma 



7r>(/3) A a 
M 



Hence (3). Conversely if ([3]) holds for a = aj{a,(^)da^ A dq^---"^ then from (Al) 

/ /3A i a= 7r^{f3)Aa = Y] / /3 A c/aM ^^ttlc^,?) 
Jn Jnj^ Jm J Jn J<;&Ma 

Since this is true for all (3 then Q holds. 

If a does not contain the factor ^^• ■^ i.e. a = aiK{<y,(i)da^ A dq^ where K ^ {1, 
then the right hand side of ([S]) becomes 

/ n'}^{(3) /\a = y2 l3jaiK{(^, ^)da-^ A da^ A dq^ = 

Jm j Jm 

and the right hand side of Q becomes 

J J ieMa 

Thus by linearity ^ and Q are equivalent for all a. 
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Lemma 3. Verification of equation 



d 4. a 



da 



Proof. Let deg(a) = p + r, deg(/9) = n—p—1 and dN and dM be the boundaries of and 
M . Since a ^ ON one may choose /3 to have support away from ON thus 



f /3 A ( i a) = 



and since a has support away from dN then 



SAT 



7rj^/3 A a = 



It follows that 



I /3a( A da) = n^{(3) A da 



N 

I d{-K'^{/3) Aa) + (-l)"-^' / d7r>(/3) A a 

-lY-r>-^ [ 7r>(/3) A « + (-1)'^-^ / n^^m A a 
JdM JM 



d/3A\ 4^ a 



j d(^(3A(^j. ay^+ j (3 Ad 
I 13 a{^1 a^ ^ I 13 Ad 



a 

a 



(3 Adi A a 
N ^J■T^r 



Lemma 4. Proof of 



-kxa = -k A a 
Px Jpx 



□ 



(A2) 



Proof. The only non-trivial a G rA{Mx x My) in (A2) can be written a = ajdx^ A dy^^"^^. 
Then 



i *x (a/rfx^ A 6/?/°^^^) = A ai{-kdx^) A dy^^'^^ = i<dx^ I ajdy^^"^^ = -k I aidx^ A dy^^"^^ 

Jpx Jpx J Mx Jpx 



□ 
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Lemma 5. H is causal on M-J if and only if 



• Z is causal on M^, 

• {dYx)\ixy) = for all {x,y) G x My such that y ^ J~{x) and 

(A3) 

• '■Sm I (^'V) ~ ^ f'^^ i^^y) ^ ^ ^Afy such that y ^ J^(a^), where 
t^My • ^ ^ ^ natural embedding. 

Proof. If isMy : ^ is the natural embedding then = ^l^^ 

and 

(xArfy(p*A))|(^,^) 



dY{tlSxAp*yA\y)\ - / i^(rfyXAp^A|,)| (A4) 

ye^MY\J-{x) ^ JyeJ^MYnJ'ix) 

^?(^^yXAp^A|,)| - /" ^^b\dYX^PYA\y)\. . 

y£M+\J-{x) ^ JyeM+nJ-{x) ^ 



First one argues that (A3) imphes that 11 is causal on Mj^. Given x G Mj^ and Fi, F2 G 
TA^M^ such that Fi\y = Fily = for ?/ G J'iy), set F = Fi - F2 so that F = on J~{x). 
Since My is topologically trivial F is exact, F = (iA, and hence dA = on J~{x). Then 
since J~{x) is topologically trivial there exists / G FA^ My such that A = df o\i J~{x). 
Thus one can choose a gauge A = A — df so that A = on J~ (x). Given <^ such that C,\y = Q 
for y G J^{x) fl Sjv/y then ^'[^jlj;. = since Z is causal. Thus from (A4) 

mc.u{x)= I ^^:^{x\^PYmy)) 



^i?4,,y (X A I + / z?,)^^^^ {x A I , ^ 



(x,y) ' / ab ^My VA l^Y \y/\(x,y) 

y&T.MY\J-{x) JyeJ^MynJ-ix) 

^^t{dYX^PYA\y)\ - [ tif{dYX^PYA\y)\,.=0 



yeM^\J-(x) JyeM+nJ-ix) 
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since l^^^ x\ix,y) = for ?/ G T,My\J {^)^ = for y E J (x) and (IyX = for ?/ G 
M+\J-(x). 

Conversely if 11 is causal on M-J then setting F = in (pi) shows that Z must be causal 



on M^. Then setting C = then for all A such that A = on J (x) (A4) yields 
= U[FXU{x) 



yeJ^MyV'ix) ^ JyeM+\J-{x) 

(A5) 

The 4-dimensional domain My\J~(x) denotes points outside the backward lightcone of 
X, while the 3-dimensional domain T,My\J~{x) denotes the points on that are not 

causally connected to x. Choosing such an A to have support about a small neighbourhood 



of y G My\J {x)\Zmy results in the first term of (A5) being zero and thus (rfy x) I (x,?/) = 0. 
Likewise setting A to have support about a small neighbourhood of ?/ G T,My\J^ (x) implies 

4My(x)l{x-,?;) = 0. 

□ 



One can now prove lemma [T] in section III B 



Proof of lemma [7} Given a E N, with V non vanishing there exists a coordinate system 
(cr^, . . . , a") on adapted to V so that V = and the image of the curve 7o- : [To(cr), 0] — )■ 
A^ is contained in the coordinate patch. Write (5 = Pjda^ then since iyP = the sum is over 
/ G {2, . . . , n}. With distinguished write /3/(cr) = (3j{a^, a) where a = (cr^, . . . , a"). Also 
since = 0, /3|(o-i,o;) = (3i{a^,a)da^ . Likewise since ivC = one has CIo-q = Ci{'^o)d(L^- 
Solving for the integral curves of V gives 0Ar(r, a^, (x) = (r + a^, a) 

and one may write TQ{a^,a) = Tq^o) — a^, giving 

^NiO\(a\a) = 0ko(<7l,a)(C|ro(<7)) = Cl {ro{a) , O) dc^ 



Thus 



( / f3i{(r^ + r,a)dT + (i{To{cT_),cr))dcr_^ 

^ 7r=ro(a)-(Tl ^ 
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Hence iyC = and one may write C\{a^,a) = 0{'^^j(L)d(L^ ■ Now 

U(^\(L) = / Pi{(T^ + r,a)dT + Cj{To{a),a) 

J T=T()(cr) — a^ 



T=TQ (ct) 



where r' = t + and 



So- J- 

9^ 



T=To(q) 



Since = on E 



N 



^\{o,a) = ^iiTo{0,a),a)da^ = ^iiTo{a), a)da^ = Ciiroia), a)da^ = C|(o,a) 



i-e. Cls^ = C- 



□ 



Lemma 6. Proof that ( 66p7 ) implies (61) and that ( 66|[69 ) implies (61) 



Proof. First (67) is equivalent to (69) since given 7 G FA^My one has i{y,u)l = u'^labQ^'^ ^ 



and hence VFL<.1(^) = ^V(j,,„)(^^) = ^u'^'jabg'"' £7 ■ From (|67|) it follows that 

x'-'Ap*y^=l^^^xdy^'A^'^:l,^^'^-^*{dTAw'^^^ Ap^7 
,w 2 



Jcl 2 



*x St^,^^'-'^{dT Awp'^ig-'^u^'^^elf)) Ap'y^ 



m 



i.e. (69). That (69) implies (67) follows since the above argument is true for all 7. 



To prove (61) note that the domains Afx and Afy'^ are related via the diffeomorphism 
TL-i :^fp^^f^\ TL"i(r,y,«) = ( - r, ^^^^^^(r)) (A6) 
Thus T^"^*{dT) = —dr and setting {x,v) = C'[^^„)(t) with r > yields 
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so that = o TL-i and thus w^^* = T^"^* o (/.L-i*. Now 



hence 



From (63) 



and from (A6) 



Hence px ° ^ = t^x ° '^x ° ^ ^^"^ 

o = TL-l'^ O t^7^1* O TT^ (A8) 

From the definition of S one has 

I S-f= I (A9) 

for any 7 G rAS(Mx x My). 

Since ^''^"^ : P — )► P' is a diffeomorphism then 

^L-i*^= / ^ (AlO) 
for any 7 G rA'^(P'). Likewise since T^"^ : A/'y^ — ?■ Afx is a diffeomorphism 



TL"i*7= / 7 (All) 



X , 



for any 7 G rA'^(A/'j 

For convenience set a\'^ = dr f\ 4''^"^* {iyi/\ci\(^p^-)OQ^') G VK^Nx ■ For fixed x assume that Fi 
has support in V^. Then one can choose /3 G FA^Mx so that p\[5 /XpyFi has support inside 



"D. Thus from (A7) 



supp(px(*/3) A = supp(p^/3 A x^"^ A pyFi) C P (A12) 
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Now 



A/x J PX 



Mx Jpx 



Mx 

JA/jf Jpx 



Mjf X My 



V 



7rx(*/3) Aa 



7ri(^/3) A i a 

X 



Mx 



>/3)A ^ A C 

' TTX 



/3 A 7^ 



Mx 



a 



TTX •J'^x 



Summing over [a] gives 



j /3A I x^PyFi = J2i'"' [ P^^l I 

JMx Jpx J Mx J T^x J "tj^ 



from (A7) 



from (A2) 



from (A9) 



from (|3]) 



from (A12) 



from (64) 



from (AlO) 



since V C M^^ 



f TL"i*w^^V^(^/3) ATL"i*aL"i from dASl) 



from (All) 



from ^ 
from ^ 



a 



Since this is true for all (3 with support in a neighbourhood of x then (61) holds at x. □ 



Lemma 7. The derivation of (72) /rom (67). 



Proof. The derivation of (72) from (67) follows by first writing the Liouville vector field (36) 



as 



^o"' = W-^-a + ^^7^ Where H'^ = -V\fU^uf + ^F.^fg'^u 



d 
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Then setting f^'^\y,u) = f^{y,u) + fl"\y,u) it follows from (31) that 

/\detg\ 



\ Uq / 



where Ya = i_a_dy^^'^^ and F = dy^^"^^. Consequently 



Un V 



.Lai |det^| 

and 

-dr A = f^^^-^^g^^vl'e^^, C^dr ^Y^ ^ du"'' + H^dr AY Adu"" 

Under the maps zuy and ip^"^* one has 

w^'^\dy'') = dy" , w^y^^du^") = du^" 

and 

^ y y J ^ ' ^ ^ ^ dy'' ^ ' dx^ dy^ 

So using the projector S* given in (70) yields 



Mo '''''' V 2 Sy'^ Sx'^ dx^ 2 dx'^ dy'' dx^ 



2 dx''- dx^ dy'^ ~^ dx'^ dx' ^ rfx 



Hence from (67) 

Jc-12 



m 



Mo V 2 (9?/'= ^x"^ 2 (9?/^ 



2 dx'^ dx^ dy^ dx'^ dx 
/ i^ildetf?! ^ /'u'^ 9r du" du'^ dr du" du" 

2 ax" ax*= oy'^ ox"- ox^ / ^ 



m 



^° 2^/0 ^ ^M-eab/^e 1^ 2 a?/'^ dx'i dx- 2 dx^ dy- dx- 



— 7—7—- h H'' — ^—]dxhi A dyi^ 

2 ax'^ ax"^ ay'^ ax'* ax'^ / 



□ 
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